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1. Introduction 



One of the most interesting problems facing the community in relativistic heavy ion collisions 
is to determine at what time the matter created can be described using hydrodynamics. In this 
context, at RHIC energies it has been found that for pr < 2 GeV, the elliptic flow of the matter 
created is described well by models which assume ideal hydrodynamic behavior starting at very 
early times T < 1 fm/c [|l|, ^, ^ ^. This is not completely understood due to the fact that the 
estimates from perturbative QCD for the thermalization time of a QGP at RHIC energies range 
from 2 — 3 fm/c [||, |6[ ^]. Moreover, recent results from conformal relativistic viscous hydrody- 
namics [||] have shown that these initial estimates for the isotropization/thermalization time of the 
plasma are not completely reliable due to poor knowledge of the proper initial conditions (CGC 
versus Glauber), details of plasma hadronization such as the choice of the proper freezeout time and 
the subsequent hadronic cascade, etc. Therefore, additional theoretical and experimental inputs are 
necessary to further constrain this time. It would be nice to have information about the thermaliza- 
tion time and appropiate initial conditions from other independent observables which are different 
than elliptic flow. One good candidate is high energy dileptons. Due to their large mean free path, 
lepton pairs can leave the strongly interacting medium after the collision. In this work, we examine 
the possibility to determine experimentally the thermalization time of the matter created using high 
energy dilepton yields as a function of both, mass and transverse momentum. We compute the 
expected e+e^ yields resulting from a Pb-Pb collision at LHC full beam energy, ^/s = 5.5 TeV in 
a pre-equilibrium scenario of a quark-gluon plasma with a time-dependent anisotropy caused by 
the rapid longitudinal expansion. We find that at LHC energies, there is an enhancement of the 
dilepton yields as a function of the transverse momentum when 2 < < 8 GeV if one chooses 
isotropization/thermalization time of 2 fm/c [|^. 

This manuscript is organized as follows: In Sec. |l] we calculate the dilepton production rate 
at leading order using an anisotropic phase space distribution. In Sec. § we review the calculation 
of the dilepton yields from an anisotropic plasma. In Sec. |3| we discuss the main aspects of the 
interpolating model from coUisionally-broadened to ideal hydrodynamical expansion. In Sec. ^we 
present the expected medium dilepton yields for LHC energies. Finally, we present our conclusions 
and give an outlook in the Sec. ^ 



2. Dilepton rate from kinetic theory 

From relativistic kinetic theory, the dilepton production rate JA'^' ' /d'^Xd^P = dR' ' /d^P at 
leading order in the electromagnetic coupling, a, is calculated as follows: [^, |l^, 11]: 

where fq q is the phase space distribution function of the medium quarks (anti-quarks), v^q is the 
relative velocity between quark and anti-quark and a^^' is the total cross section. Since we will be 
considering high-energy dilepton pairs with center-of-mass energies much greater than the dilepton 
mass we can ignore the finite lepton mass corrections and use simply a^^' = 47ra^/3M^. In 
addition, we assume that the distribution function of quarks and anti-quarks is the same, fq = fq. 
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We will follow the ansatz by Romatschke and Strickland for the distribution function that 



describes the pre-equilibrated stage of the plasma [|12|]. The ansatz consists in consider that any 
anisotropic distribution function is obtained from an arbitrary isotropic phase space distribution by 
squeezing {t, > 0) or stretching {t, < 0) the isotropic distribution function: 



(2.2) 



where phard is the hard momentum scale, fi is the direction of the anisotropy and ^ > is a parameter 
that reflects the strength and type of anisotropy. In the isotropic case, phard is identified with the 
temperature of the system and ^ =0.ln this work we will study the case when the direction of the 
anisotropy is along the longitudinal (beamline) direction, i.e., n = e^. This configuration is relevant 
for earliest stages of the collisions between heavy nuclei. 
From Eqs. (2.2) and (2.1) we obtain 



dR 



5a' 



J4p ~ ISn^ .L/^^"^^P'\fa /5(^^pf(l + i^(T)cOs2 0pJ,;^hard(T,T7: 

X (^^ (E - pi)2 + ^ (t) (pi COS - PCOS 0p)2,Phard(T, T] ) j , (2.3) 



with 



^ = 4P';?f sin' Opsin' 0p, -(2pi(£- P cos 0p cos 



pi- 



-M 



2\2 



M' 



2(£-Pcos(0p±0p,)) 



Notice that the dilepton rate does not tell us anything about the space-time evolution of the system 
and, therefore, it is not enough to make a phenomenological prediction for the expected dilep- 
ton yields. In order to make contact with experiments, it is necessary to include the space-time 
dependence of /?hard and ^ and then integrate over the space-time volume 



dN 



nr 



dM^dy 

dN^^^' 
d^prdy 



kRj 



d^PT 



dM^ 



To v'-c 



dR 



. d^P 
dR^ 
d'^P 



-Tdrdrj 



-Tdzdrj 



(2.4a) 
(2.4b) 



where Rt = 1.2A^/^ fm is the radius of the nucleus in the transverse plane. These expressions are 
evaluated in the center-of-mass (CM) frame while the dilepton production rate is calculated for the 
local rest frame (LR) of the emitting region. Then, the dilepton pair energy has to be understood 
as Elr = pj cosh(3' — T]) in the differential dilepton rate dRann/d^P. Substituting Eq. (23) into 



Eqs. ( |2.4t ), we obtain the dilepton spectrum including the effect of a time-dependent momentum 



anisotropy. In Eqs. ( |2.4| ) it is assumed that there is only longitudinal expansion of the system. 
This assumption is well justified since corrections to high energy dileptons coming from transverse 
expansion or mixed/hadronic phase do not play an important role in the studied kinematic regime 



[14] 



'Details of the calculation are presented in Ref. [p^. 
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Note that we have not included the next-to-leading order (NLO) corrections to the dilepton 
rate due to the complexity of these contributions for finite t, . These affect dilepton production for 
isotropic systems for < 1 [[11, |T6[ |T7|, [Isl]. In the regions of phase space where there are large 
NLO corrections, we will apply ^-factors to our results as indicated. 



3. Space-time model 

Previous phenomenological studies of high energy dileptons have assumed that the value of the 
isotropization time, Tjso, is of the same order as the parton formation time Tq. However, the physical 
mechanisms which could make such fast isotropization feasible are not completely understood. 
Recently, a phenomenological model has been proposed where the isotropization time could be 



larger as a consequence of a previous evolution in a pre-equilibrium dynamics [[131, |14|]. Before 
going into the details of the model, we remind the reader of some general relations. 

First, the plasma anisotropy parameter t, is related with the average longitudinal momentum 
{pi) and transverse momentum {pj) of the hard particles through the expression: 

Note that this relation is completely general and can be applied in all cases. It is possible to obtain 



two important limits from Eq. (3J_). When we have that {pj) = 2{p\), then § = 0. This is the case 
when the system is isotropic in momentum-space. Another possibility is when the partons expand 
freely along the longitudinal axis, i.e.. Id free streaming expansion. Using the free streaming 
distribution function, it is possible to show that the transverse and longitudinal momentum scales 



as [|, |g,|l9p: 

{p'r)t,.oc2T^ , (3.2a) 



2 1^2 _ 



{Ptk.. - 7? -I . (3.2b) 



Inserting these expressions into Eq. jpA]), one obtains i^y.^. (t) = T /Tq — 1. The free streaming re 



suit for ^ is the upper-bound on possible momentum-space anisotropics developed during Id expan- 



sion by causality. When the system has different kind of interactions, Eq. (|3JJ) for the anisotropy 
parameter will scale differently as discussed below. 

Second, for a given anisotropic phase space distribution of the form specified in Eq. (p^), the 
local energy density can be factorized as: 

^(/^h^d, = J -0^ P Ao(v'p2 + §(p-A)2,;^hard) , (3.3) 

= MPhard)^{^), 

where Sq is the initial local energy density deposited in the medium at Tq and 

^/t-x 1 / 1 arctanA/?\ 

In this work, we will not study explicitly the possibility of Id free streaming expansion since, in 
reality, this is a rather extreme assumption which requires that the partons do not interact at all [pj|]. 
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3.1 Momentum-space broadening and plasma instabilities effect 

The ratio between the average longitudinal and transverse momentum needed to compute 
using Eq. ( ^ ) is modified from the free streaming case if collisions between the partons are taken 
into account. In general, it is a difficult task to obtain the exact form of collisional kernel of the 
Boltzmann equation plus mean field interactions (Vlasov term). As a first approach, one can start 
considering the effect of elastic collisions in the broadening of the longitudinal momentum of the 

particles. This was the approach in the original version of the bottom-up scenario [S]. In the first 

3/2 

Stage of this scenario, 1 ^ QsT <^OCs , initial hard gluons have typical momentum of order and 
occupation number of order 1 /a^ . Due to the fact that the system is expanding at the speed of light 
in the longitudinal direction, the density of hard gluons decreases with time, Ng ~ / {cCsQ^x) . 
If there were no interactions this expansion would be equivalent to 1+1 free streaming and the 
longitudinal momentum pi would scale like 1 /t. However, once elastic 2^2 collisions of hard 
gluons are taken into account, the ratio between the longitudinal momentum pi and the typical 
transverse momentum of a hard particle pr decreases as : 



{P 



2\ 



(Pi) 



{QsT)-^'^ . (3.5) 



Assuming, as before, isotropy at the formation time, Tq = Qs ^ this implies that for a coUisionally- 
broadened expansion, (^(t) = (t/Tq)^^^ — 1- 



Note that in order to obtain the relation given by Eq. ( |3.5[ ) it is implicitly assumed that that 
the elastic cross-section is screened at long distances by an isotropic real- valued Debye mass |^]. 
But this is not the case of an anisotropic plasma, since the Debye mass can become complex due to 
the chromo-Weibel instability |[l^, |l^, 20, 21], 22]. Such negative eigenvalues indicate instabilities, 
which result in exponential growth of chromo-electric and magnetic fields, and B", respectively. 
These fields deflect the particles and how much deflection occurs will depend on the amplitude and 
domain size of the induced chromofields. Currently, the precise parametric relation between the 
plasma anisotropy parameter and the amplitude and domain size of the chromofields is not known 
from first principles. 

If one would like to include the momentum-space broadening and the impact of the plasma 
instabilities, this can be achieved by considering the temporal dependence of t, (t) as: 

^(T,5)=(|y-1. (3.6) 

The exponent 5 contains the physical information about the particular type of momentum-space 
broadening which occurs due to plasma interactions. Two limiting cases for this exponent are the 
ideal hydrodynamic and free streaming expansion. In the l-i-l hydrodynamical limit, 5 = and then 
1^ ^ 0. For 5 = 2, one recovers the l-i-l dimensional free streaming case, t, —>■ i^f s. = (t/tq)^ — 1- 
For < 5 < 2, one obtains the proper-time dependence of the energy density and temperature 
by substituting ( ^^ ) into the general expression for the factorized energy density ( pj| ) to obtain 

<f(T,5)=^0(;^hard)=^^(^(T,5)). 

Different values of 5 arise dynamically from the different processes contributing to parton 



isotropization. For example, elastic collisional-broadening results in Eq. (3.5) and hence 5 = 2/3 
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Recently, some autiiors iiave considered the values of 5 resulting from processes associated with the 



chromo-Weibel instability presented at the earliest times after the initial nuclear impact [ )23i gj, |25|] : 



where 



Ref.[231 



1 Ref. [24] 



(3.7) 



(3.8) 



2 Nielsen-Olesen limit, Ref. [25] . 



Summarizing, the coefficient 8 in Eq. (3.6) takes on the following values 



2 Free streaming expansion , 
2/3 CoUisional-Broadening, Ref.[ 



1/2 Ref.[23] , 



1/4 Ref.[24] , 



(3.9) 



1/6 Nielsen-Olesen hmit, Ref.[g5|] , 
Hydrodynamic expansion . 

In Fig. [T] we sketch the time-dependence of the plasma anisotropy parameter indicating the 
time scales at which the various processes become important. At times shorter than the mean 
time between successive elastic scatterings, Tmfp. the system will undergo l-i-l dimensional free 
streaming with 8 = 2. For times long compared to Tmfp but short compared to Tinstabiiity the plasma 
anisotropy will grow with the collisionally -broadened exponent of 8 = 2/3. Here ^instability is the 
time at which instability-induced soft gauge fields begin to influence the hard-particles' motion. 
When Tinstabiiity < T < Tiso the plasma anisotropy grows with the slower exponent of 5 = 1 /6 . . . 1 /2 
due to the bending of particle trajectories in the induced soft-field background. At times large 
compared to Tinstabiiity inclastic processes are expected to drive the system back to isotropy [j5|]. We 



note here that for small and realistic couplings it has been shown [26] that one cannot ignore the 
effect of coUisional-broadening of the distribution functions and that this may completely eliminate 
unstable modes from the spectrum. 

From this sketch, it may be possible to try to construct a detailed model which includes all 
of the various time scales and study the dependence of the process under consideration on each. 
However, there are theoretical uncertainties in each of these time scales and their dependences 
on experimental conditions. Because of this, we will construct a phenomenological model which 
smoothly interpolates the coefficient 8 from the Id coUisionally-broadened expansion to Id hydro- 
dynamical expansion, i.e., 2/3 > 5 > 0. 

In the model we introduce a transition width, y^^, which indicates the smoothness of the tran- 
sition from the initial value of 5 = 2/3 to 5 = at T~ Tjso. Note that by using such a smooth in- 
terpolation one can achieve a reasonable phenomenological description of the transition from non- 
equilibrium to equilibrium dynamics which should hopefully capture the essence of the physics. 
The coUisionally-broadened interpolating model provides us a realistic estimate of the effect of 
plasma anisotropics. 
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Figure 1: Sketch of the time dependence of the plasma anisotropy indicating the various time scales and 
processes taking place. Here Tmfp is the mean time between elastic collisions (mean-free time) and Tinstabiiity 
is the time at which plasma-instability induced soft modes have grown large enough to affect hard particle 
dynamics. 



3.2 Interpolating model for coUisionally broadened expansion 



In order to construct an interpolating model between coUisionally-broadened and hydrody- 
namical expansion, we introduce the smeared step function: 



A(T,Tiso,7) = ^ (tanh 



+ 1 



(3.10) 



where y ' sets the width of the transition between non-equilibrium and hydrodynamical evolution 
in units of Tiso. In the limit when T <^ Tjso, we have A ^ and when T » Tjso we have A ^ 1. Phys- 
ically, the energy density S should be continuous as we change from the initial non-equilibrium 
value of 8 to the final isotropic 5=0 value appropriate for ideal hydrodynamic expansion. Once 
the energy density is specified, this gives us the time dependence of the hard momentum scale. We 
find that for general 8 this can be accomplished with the following model for fixed final multiplicity 
Pll: 



^K5) 

/?hard(T,T]) 



(3.11a) 
(3.11b) 
(3.11c) 
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with ) defined in Eq. ( p] ) and for fixed final multiplicity we liave: 

^(T) = [i^((TisoAo)^-l 

^ lini ^(t) = 



3A(T)/4 ^Tiso\ 1-5(1-^ 



(Tiso/^o) -1 



To 



(3.12a) 



(3.12b) 




and 5 = 2/3 for the case of Id collisionally broadened expansion interpolating to Id ideal hydro- 
dynamic expansion. 
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Figure 2: Temporal evolution using our fixed final multiplicity interpolating models for the energy density 
(left column), hard momentum scale (middle column), and anisotropy parameter (right column) for four 
different isotropization times Tiso G {1,4,6, 18} Tq. To convert to physical scales use Tq ^ 0.3 fm/c for RHIC 
and To -0.1 fm/c for LHC. 

In Fig. ^ the temporal evolution of the energy density, hard momentum scale and the anisotropy 



parameter ^{t) is plotted using Eq. (3.11a). As can be seen from the figure for fixed final multi- 
plicity at T = To, the initial value of Phard is reduced for finite values of Tiso but once the system 
looks isotropic in momentum-space around T > Tiso, Phard is the same, independent of its initial 



value [13]. 



4. Results 

We calculate predicted e^e^ yields as a function of invariant mass and transverse momentum 
along with predicted yields from other sources for LHC energies. For comparison with previous 
works we take To = 0.088 fm/c, Tq = 845 MeV, = 160 MeV, and Rt = 7.1 fm Here we 
assume that when the system reaches Tc all medium emission stops. 

In Fig. m we show our predicted dilepton mass spectrum for LHC energies using the time 
dependence of the energy density, the hard momentum scale and the anisotropy parameter given by 
Eqns. (3.11) with 5 = 2/3. As can be seen from Fig. |3| there is a difference of the medium dilepton 
yield when varying the assumed plasma isotropization time from 0.088 fm/c to 2 fm/c. However, 
the prediction is up to one order of magnitude below the other contributions of dilepton yields 
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Figure 3: CoUisionally-broadened interpolating model dilepton yields as a function of invariant mass in 
central Pb+Pb collisions at the LHC, with a cut pj > 8 GeV and and rapidity 3'=0. For medium dileptons 
we use 7=2 and Tiso = 0.088 fm/c for LHC energies and fixed final multiplicity. A /T-factor of 1.5 was 
applied to account for NLO corrections. Dilepton yields from Drell Yan, Heavy Quarks, Jet-Thermal and 
Jet-Fragmentation were obtained from Ref. 



(Drell- Yan, jet-thermal, and jet-fragmentation). This coupled with the large background coming 
from semileptonic heavy quarks decays would make it extremely difficult for experimentalists to 
extract a clean medium dilepton signal from the invariant mass spectrum. For this reason it does 
not look very promising to determine plasma initial conditions from the dilepton invariant mass 
spectrum. Nevertheless, the situation looks better for the dilepton spectrum as a function of the 
transverse momentum as it is shown in Fig. 0. In this spectrum, there is an enhancemment of 
medium dilepton yield for 2 < pr < S GeV compared with the other sources of dileptons (Drell 
Yann and Jet conversion). Therefore, this observable offers the oportunity to have a cleanest way 
to determine plasma initial conditions. 

In order to quantify the effect of time-dependent pre-equilibrium emissions we define the 
"dilepton enhancement", 0(Tiso), as the ratio of the dilepton yield obtained with an isotropization 
time of Tiso to that obtained from an instantaneously thermalized plasma undergoing only l-i-l 
hydrodynamical expansion, ie. Tiso = '^o- 



0(Tiso) = 




(4.1) 



This ratio measures how large the effect of early-time momentum anisotropics are on medium 
dilepton production. In the case of instantaneous isotropization, <I>(Tiso) is unity, and for Tiso > 
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Figure 4: Collisionally -broadened interpolating model dilepton yields as a function of transverse momentum 
in central Pb+Pb collisions at the LHC, with a cut 0.5 < M < 1 GeV and rapidity y=0. For medium dileptons 
we use 7=2 and Tiso = 0.088 fm/c for LHC energies and fixed final multiplicity. A A'-factor of 6 was applied 
to account for NLO corrections. Dilepton yields from Drell Yan, Jet- Thermal and Jet-Fragmentation were 
obtained from Ref. [111. 



To any deviation from unity indicates a modification of medium dilepton production due to pre- 
equilibrium emissions. 

In Fig. ^ we show the dilepton enhancement, , as function of transverse momentum for Tiso = 
2 fm/c at LHC energies. The invariant mass cut is the same as in Fig. ^ (0.5 < M < \ GeV). As can 
be seen from Fig. ^ there is a rapid increase in 1 and 4 GeV at LHC energies. Moreover, from 
this figure both sharp and smooth transitions from early-time collisionally-broadened expansion 
to ideal hydrodynamic expansion result in a 30-50% at LHC energies. The effect of reducing Tiso 
is to shift the peak in to larger pT while at the same time reducing the overall amplitude of 
the peak. Therefore, in order to see the difference between an instantaneously thermalized QGP 
with Tiso = ^0 and one with a later thermalization time requires determining the medium dilepton 
spectra between 2<pt <^ GeV at LHC with high precision so that one could measure the less than 
50% variation resulting from pre-equilibrium emissions. Finally, we point out that it is possible to 
take other cuts (invariant mass and/or transverse momentum). This could be coupled with fits to 
experimental data, allowing one to fix Tiso and y via a "multiresolution" analysis. 



5. Conclusions 



In this work we have introduced a phenomenological model that takes into account early-time 
momentum-space anisotropics in the rapidity dependence of high-energy dilepton production. To 
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Figure 5: Dilepton enhancement, as defined in Eq. (4.1) resulting from our coUisionally-broadened in- 
terpolating model with fixed final multiplicity and 7=2. The invariant mass cut used was 0.5 < M < 1 
GeV, rapidity ^=0 and To = 0.088 fm/c for LHC energies. Lines show expected pre-equilibrium dilepton 
enhancements for different values of the assumed plasma isotropization time, Tiso- 



do this we have modeled the temporal evolution of the plasma anisotropy parameter and the hard 
momentum scale phard- Using these models we integrated the leading order rate for dilepton pro- 
duction in an anisotropic plasma over our modeled space-time evolution. Based on our numerical 
results for the variation of dilepton yields with the assumed values of Tjso, we find that the best 
opportunity to determine information about the plasma isotropization time is by analyzing the high 
transverse momentum (2 < /jj < 8 GeV at LHC) dilepton spectra using relatively low pair invari- 
ant mass cuts (M < 2 GeV). Based on these pj spectra we introduced the "dilepton enhancement" 
factor (^>(Tiso) which measures the ratio of yields obtained from a plasma which isotropizes at Tjso to 
one which isotropizes at the formation time, Tq. We find that there is a 30-50% enhancement in the 
high-transverse momentum dileptons at LHC energies when one assumes an isotropization time of 
Tiso = 2 fm/c. The amplitude of the enhancement and position of the peak in the enhancement func- 
tion, 0, varies with the assumed value of Tiso which, given sufficiently precise data, would provide 
a way to determine the plasma isotropization time experimentally. Another interesting possibility 
to constrain the value to Tiso using dileptons is studying the forward region, where the effect of 
early-time anisotropics is expected to be maximal |27]. 

As an extension of the interpolating model presented in this work, one can consider the case 
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where, instead of having at late times 1+1 ideal hydrodynamical expansion, the system is described 
during that stage by 1+1 viscous hydrodynamical expansion. This is possible since one can relate 
the anisotropy parameter, , with the anisotropy in momentum-space of the fluid generated by the 
difference of the pressures along the longitudinal and transverse plane due to viscous corrections 

ii. 

An uncertainty of our treatment comes from our implicit assumption of chemical equilibrium. 
If the system is not in chemical equilibrium (too many gluons and/or too few quarks) early time 
quark chemical potentials, or fugacities, will affect the production of lepton pairs [|l^, [IT|]. How- 
ever, to leading order, the quark and gluon fugacities will cancel between numerator and denomina- 
tor in the dilepton suppresion factor, 0(Tiso) [11 1. We, therefore, expect that to good approximation 
one can factorize the effects of momentum space anisotropics and chemical non-equilibrium. 

We note in closing that within this model it is possible to assess the phenomenological con- 
sequences of momentum-space anisotropics on other observables which are sensitive to early-time 
stages of the QGP, e.g. photon production, heavy-quark transport, jet-medium induced electromag- 
netic and gluonic radiation, etc. 
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